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r S^ . Abstract 
-i— > . 

We study statistical models, specifically transfer matrices corresponding to a mul- 
tiparameter hierarchy of braid matrices of (2n) 2 x (2n) 2 dimensions with 2n 2 free 
parameters (n = 1, 2, 3, . . .). The simplest, 4x4 case is treated in detail. Powerful 
recursion relations are constructed giving the dependence on the spectral parameter 
■ 9 of the eigenvalues of the transfer matrix explicitly at each level of coproduct se- 

■"sj- ! quence. A brief study of higher dimensional cases (n > 2) is presented pointing out 

features of particular interest. Spin chain Hamiltonians are also briefly presented for 
the hierarchy. In a long final section basic results are recapitulated with systematic 
analysis of their contents. Our eight vertex 4x4 case is compared to standard six 
5o ' vertex and eight vertex models. 
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(N 

o 



1 Introduction 

In a previous paper p] statistical models were presented starting from a class of multipa- 
rameter braid matrices of odd dimensions (Ref. 1 cites previous sources). This class of 
N 2 x N 2 braid matrices for iV = 2n - 1 (n = 2, 3, . . .) depends on \ {N + 3) (N - 1) = 

2 (n 2 — 1) parameters and has a basis of a nested sequence of projectors introduced before 
[2]. Recently this class has been extended to include even dimensions [3]. The even di- 
mensional (2n) 2 x (2n) 2 braid matrices depend, apart from the spectral parameter 9, on 
2ra 2 free parameters. One such parameter can, as usual, be absorbed by a suitable choice 
of an overall normalization factor. Such a quadratic increase of free parameters with n 
(even or odd) is the most salient, unique feature of our models with their nested sequence 
basis. 
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We start with a detailed study of the 4x4 case. This is our new eight vertex model. 
It should be compared to exotic eight vertex model based on the S03 bialgebra [I]. There 
negative Boltzmann weights were involved and were commented upon. Here, for suitable 
choices of parameters, one has all weights positive or zero for 9 > or for an alternative 
choice for 9 < 0. This 4x4 case has another major distinguishing feature. It finds its 
place as the first step in an explicitly constructed multiparameter hierarchy. This is to be 
contrasted with the multistate generalization of six vertex model (Ref. 5 and the original 
sources cited there) in which the parametrization remains restricted to the six vertex level 
as the dimensions increases. We will present some remarkable general features conserved 
in the entire hierarchy. But here our study will be brief - a beginning. We hope to explore 
further elsewhere. We will also briefly present spin chain Hamiltonians for our class of 
solutions. Some aspects are better discussed after the constructions are presented. Such 
comments are reserved for the end. 



2 Even dimensional multiparameter braid matrices 



We briefly recapitulate the construction of Ref. 3, closely related to that for odd di- 
mensions [U [2]. The notation of Ref. 3 are maintained. Define the nested sequence of 
projectors 

1 



P 



{(ii) ® (jj) + (ii) <g> (jj) + e [(ii) ® (jj) + (ii) <g> (jj))} , 



(2.1) 



where i, j 6 {1, • • • , n}, e = ±, i = In + 1 — i, j = 2n + 1 — j and (ij) denotes the matrix 
with an unique non-zero element 1 on row i and column j . They provide a complete basis 
satisfying 



n 



! (2n) 2 x(2n) 2 , 



(2.2) 



where P.-) is obtained by changing j with j in P^f '. The (2n) 2 x (2n) 2 braid matrices 



(e) 



e=± i,j=l 



p( e ) , p(f) 



satisfy 



R12 (0) R 23 (9 + 9')R 
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R23 {9')Ri2 {9 + 9') R 23 (9) 



(2.3) 



(2.4) 



where R\2 = R®I and R23 = I®R. Here, as for odd dimensions [Hl2] the free parameters 
enter as exponents (^e m ^ and the basic constraint [H [2] = m^) has already been 
incorporated in (2.3). This leaves 2n 2 free parameters. The simplest examples are: 
• jV = 2n = 2: 



R{9) 
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(2.5) 
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where 



• N = 2n = 4: 
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/ 
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22 



22 



A22 — A- 



22 



and 



2 



6± 
d ± 



1 Cem^e ± e m^e 



5( 



e m 22 )fl -t- g^L' 



a± = 

c± = ±(e m & )9 ±e m k )e '), 

We will often write 

Dn = (a+, K, K, a +)diag. ' A ^ = ( a -' a -)anti-diag. 
and so on. Generalization for n > 2 is entirely straightforward. 



(2.6) 



(2.7) 
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(2.8) 



(2.9) 



(2.10) 



3 i?TT relations and transfer matrix (n = 1) 

To simplify computations, taking out an overall factor a + we write R (9) of (2.5) as 



i?(x) 



1 








X 
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X 








X 


1 





X 








1 



(3.1) 



where 



X = ~^h> ~» = tanh ^ ' 



(3.2) 
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with 2fi 



m^ 1 — m![i ■ Note that 



jmft) > mfc\ 9 > 0} =^x > and {m^ < m^ , < o} x > 0. (3.3) 

We consider these two domains separately, assuring nonnegative, real Boltzmann weights. 
For purely imaginary parameters im^ (m^ real) 



x = itan (fi9) , 

with a normalization factor cos fi9, one obtains a unitary braid matrix 

J R(x) + i?(x) = /. 



(3.4) 
(3.5) 



This has been pointed out in sec. 3 of Ref. 3 citing Ref. 6. We continue to study here 
real matrices. 

The basic _RTT equation can now be written as 



R (x") (T (x) ® T (x')) = (T (x') <S> T (x)) R (x") , 



(3.6) 



where 



x = tanh (/jl9) , x' = tanh {jJ,9') 



x — x 



X 



1 - XX 

We denote T = T (x), V = T (x'), R" = R (x") and the 4 blocks of T as 



tanh {9-9'). (3.7) 



T = 



Tu 


Ti 2 
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T 21 


T22 




C 





Defining the 2x2 matrices 



1 
1 



K 



1 

1 



(3.8) 



(3.9) 



(3.6) is now (with A = A (x), A' = A (x'), etc.) 
(I 01 + x."K <g> X) 









A' 


5' 




A' 


B' 
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c 






C" 
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D' 
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D 



(I®I + x."K®K). (3.10) 



A detailed study of (3.10) is given in Appendix. Recursion relations will be extracted from 
it and implemented to construct the spectrum of the eigenvalues of the transfer matrix in 
sec. 5. The transfer matrix (A + D) r = A r + D r is obtained by starting with fundamental 
2x2 blocks (for n — 1) and constructing 2 r x 2 r blocks using standard prescriptions 
(coproduct rules - see sec. 5). The starting point, the 2x2 blocks for r = 1, is provided 
by the Yang-Baxter matrix corresponding to (3.1), namely 



R (x) = Pi? (x) 



1 x 
0x10 
1x0 
x 1 



A l B l 
Ci D 1 



(3.11) 
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where P is the standard 4x4 permutation matrix. Then one uses the coproduct rule 

( T ^) r+1 = E( T ^i®( T ^)r- ( 3 - 12 ) 

k 

Such a construction guarantees the commutativity of the transfer matrix , i.e. 

[T(x),T(x')] = (3.13) 

quite generally and reducing for our present case to 

[(A + D) r , (A + D)' r ] = (3.14) 

for all r. This is the basic ingredient of exactly solvable statistical models [7]. The 
trace and the highest eigenvalue of T (x) provide significant features of the corresponding 
models. We will obtain these quite simply and generally for our class. For n > 1 certain 
essential features will be briefly presented in sec. 6. 



4 Eigenfunctions and eigenvalues of the transfer ma- 
trix 

The next essential step is the construction of eigenstates and eigenvalues of the transfer 
matrix. In this section we display the results for r = 1,2,3,4. They provide explicit 
examples of the iterative structure of the transfer matrices T r to be derived in the next 
section. Moreover they illustrate how multiplets involving r-th roots of unity and pos- 
sible multiplicities of them combine to provide a complete basis of mutually orthogonal 
eigenstates of T r spanning the base space. Referring back to them one grasps better the 
full content of the general formalism of sec. 5. We start with the notations of sec. 3. 

• r = 1: 



T 1 = (A + J D) 1 = (l + x) 
The eigenstates are evidently 

with 



1 
1 



(A-D\ = (l-x) 



1 
-1 



(4.1) 



|1> 



and 



r = 2: 



T 1 (|l),|l» = (l + x)(|l),|l» 
Tr(T 1 ) = 2(l+x). 

T 2 = (A + D) 2 = (1 + x) 2 X (2 ,o) + (1 - x) 2 X(o, 2) , 



(4.2) 

(4.3) 
(4.4) 

(4.5) 
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where 

X( 2 ,0) = - 



1 






® 


i 





2 











1 



with I 



(i) 



1 
1 



and -£T(i) 



X( ,2) = - 

One has 

X(2,0)X(o,2) 



1 


hi) 






1 


2 




-hD 




-1 



(4.6) 



1 

1 



X(0,2)X( 2 ,0) = 0. 



(4.7) 



We have anticipated the iterative structure of sec. 5. That T 2 obtained by a straightfor- 
ward use of coproduct rules as 



1+x 2 



2x 





2x 

;i+x) 
o 



2x 

2 



;i + x 

2x 







1 + x 2 



(4.8) 



can be expressed as (4.5) on a basis satisfying (4.7) is the central lesson. Denote 



|H), 



111), 



111), 



|11). 
(4-9) 

For r > 1 the order of the indices (1, 1) indicates the structure of the tensor product. For 

^ ). One obtains, with upper or lower signs, 



example, 1 1 1 1 ) 







J)® 





T 2 (|11) ± |II)) = (1 ± x) 2 (|11) ± |II)) , T 2 (|1I) ± |I1)) = ± (1 ± x) 2 (|1I) ± |I1)) 

(4.10) 

and 

Tr(T 2 ) = 2 (1 + x) 2 . (4.11) 

The eigenfunctions of X( 2 ,o) (resp. X( 0;2 )) are annihilated by X( 0;2 ) (resp. X( 2j o)) consis- 
tently with (4.7). Finally, again anticipating sec. 5, 



(A-£>) 2 = -(l+x) 



hi) -Kq.) 
~ K (i) hi) 



(A-D) 1 + -(l-x) 



hD 



-hD 



(A + D\. 
(4.12) 

• r = 3 : From this stage onwards one can better appreciate the role of multiplets involving 
roots of unity in constructing eigenstates. One obtains by implementing coproducts in 
step (r = 2) — ► (r = 3), 

T 3 = (A + D) 3 = (1 + x) 3 X (3 ,o) + (1 + x)(l - x) 2 X (1)2) , (4.13) 

where the 8x8 matrices satisfy 



X(3,o)X( lj2 ) = Xn ^X 



(1,2)^(3,0) 







(4.14) 
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and they are x-independent. They are easy to obtain but will not, for brevity, be presented 
explicitly. The base space splits up into 4-dim. ones, closed under the action of T 3 and are 
characterized by even (odd) number of the index 1 (considering zero as even) respectively. 
Define 

Id) = |111) + |111) + |111) + |111), 
|e 2 ) = -3 |m) + 1 11 1> + |111> + 1 11 1> , 
|e 3 ) = |lll) + w|lll)+cj 2 |lll), 

|e 4 ) = 1 111> + u 2 |111> + u |111) , (4.15) 

where to = e 1 ^" (oj 3 = 1 and 1 + uj + uj 2 = 0). One has (e^e,) = 0, i ^ j, where (ej| 
denotes the transform of |ej) with conjugated coefficients. Denote also 

Id) = |111) + |111) + |111) + |111), 
|o 2 ) = -3 |111) + |111) + |TlT) + |111) , 
\o 3 ) = 1 111> + + uj 2 |111) , 

|o 4 ) = |111> + to 2 |m) + uj |111) . (4.16) 

The states {|ej) , \oj)} form a complete basis of orthogonal states. Consistently with 
(4.13), (4.14) one obtains 

T3| ei ) = (l + x) 3 | ei ), 

T 3 (|e 2 ) , |e 3 ) , |e 4 )) = (1 + x) (1 - x) 2 (|e 2 ) , u |e 3 ) , oo 2 |e 4 )) , 
T 3 |oi) = (l+x) 3 |o!>, 

T 3 (|o 2 ) , \o 3 ) , |o 4 » = (1 + x) (1 - x) 2 (|o 2 ) , oo |o 3 ) , oo 2 |o 4 » (4.17) 

and, finally, 

Tr (T 3 ) = 2 (1 + x) 3 + 2 (1 + x) (1 - x) 2 (l + u + uj 2 ) = 2 (1 + x) 3 . (4.18) 
• r = 4 : Now 

T 4 = (A + D) 4 = (1 + x) 4 X ( 4,o) + (1 + x) 2 (1 - x) 2 X (2)2) + (1 - x) 4 X (0 ,4), (4.19) 
where the 16 x 16 constant matrices satisfy 

X(4,0)X(2,2) = X( 4)0 )X(o,4) = X( 2j 2)X(o ) 4) = X( 2j2 )X( 4i0 ) = X( 0l 4)X( 4)0 ) = X(o,4)X( 2j2 ) = 

(4.20) 

as consequence of recursion relations involved in (A ± D) 3 — >■ (A ± D) 4 . The X's are 
obtained fairly easily. Now the even and odd subspaces are 8-dim. We display the 
eigenstates explicitly to illustrate a new feature. Now r = 4 is not a prime number and 4- 
th roots and square roots of unity (corresponding to r = 2 x 2) both contribute multiplets 
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involving respective coefficients (1, i, — 1, — i) and (1, —1) (For r = 6 one would have thus 
2-plets, 3-plets and 6-plets). The |e) and |o) spaces have the orthogonal bases 



|ei> = 


|iiii) + 


|im) + |iiii) + 


|1111) + |1111) + |1111) 


+ |iiii) + |iiii), 


h) = 


|iiii) - 


|im) + 1 1 Iii> - 


|im) - |i n i> + | iiii> 


-inn) + |iiii), 


h) = 


|iiii) - 


|iiii), 






M = 


|iiii) - 


|iiii) - |iiii) + 


|iiii), 




N) = 


|iiii) - 


|iiii) - |iiii) + 


|iiii), 




N) = 


|im) - 


|mi) , 






M = 


|im) - 


i 1 1 1 11> + i |1111) 


-inn), 




h) = 


|iiii) + 


i |1111) — i ]1H1> 


-inn), 


(4.21) 



and 



h) 

\°2) 

l<*> 
M 

l°6> 
\<*) 
\or) 
\os) 



1111) + |llll) + |llll) + |Illl) + |1III) + |1111) + |1111) + |1111) , 
1111) - |1111) + |1111) - |1111) + |1111) - |1111) + |1111) - |1111) , 

mi) + | ii ii> + |im) + |Im) - |mi) - |iiii> - |ml) - |im) , 
mi) - |nii) + |iin) - |iin) - |iiii> + |iiii) - |HII) + 1 1 ill) , 

1111) - i |1111) - |1111) + i |1111) , 
1111) + i |Im) - 1 1 11 1> - i 1 11 11> , 
1111) + i|llll) - | 111 1> — i lllll) , 

1111) - i |1111) - |1111) + i |1111) . (4.22) 



Define 



T 4 \a) = vf ] , T 4 \ 0i ) = \ 0i ) , i = 1, . . . , 8. 
Then one obtains, in order, 



(4.23) 



(e) 

,(°) 



(e) 
8 



(1 + x) 4 , (1 - x) 4 , (1 - x) 2 (1 + x) 2 (1, 1, -1, -1, i, -i) 

\4 /i \4 \2 / n , \2 



and 



r, ... .r, =(l + xr,-(l-xr,(l-x)^(l+x)^(l,-l,i,i,-i,-i) (4.24) 

Tr (T 4 ) = 2 (1 + x) 4 . (4.25) 



5 Relating (A, B,C,D) for all r and constructing iter- 
atively the eigenvalue spectrum 

2, to the 



We consider below exclusively the 4x4 case. Coproduct rules lead, for N 
recursion relations 



A r+ i 
B r+ i 



A, 



A r + Bi 
B r + B\ 



D r+ 1 = D 1 

C r+ i = C\ 



D r + Ci 



>B r 
C r 



(5.1) 



where 

Thus 

A r+ i 

Denote i 



A 1 



1 
x 



D 1 



x 
1 



D 



r+l 



x£L D r 



1 
1 



, K = K 2 = 



, B 

1 

1 



r+l 



x 

1 



B r xD r 
D r xIL 



Ci 



r+l 



1 
x 



XC^ 



and for p factors 1^ 



(5.2) 



(5.3) 



I, 



= K (g) K (g) ... (g) K. Starting from B x = KA U C x = KD U D x = KA X K and 
iterating one obtains 



B r = (l( r -l) <g)K)A r , C r = (l( r -l) <S> K) D r 

D r = K {r) A r K {r \ A r = K {r) D r K {r) . 



(5.4) 



Using (i(r-i) <E> if) = if ^ 1 - ) <E> J one can express (A, 5, C, Z?) r in term of any one of 
them. In particular, 

(B r ± C r ) = (J (r _i) <g> if) (A ± A.) , (5.5) 

where 

ftr-D ® if ) = ^ • • • , if ) d iag. = *M ( 5 - 6 ) 

Thus 

if 
if 
if 
if 

Thus exchanging members of successive pairs of rows [(1, 2), (3, 4), . . . , (2p — 1, 2p), . . -]( rows ) 
one obtains (B r ± C r ) from (A r ± D r ). The following recursion relations are implied (with 
the definition (5.6) of if( r )) 



(B 3 ± C 3 ) 



(A3 ± £> 3 ) . 



(5.7) 



(A + L>) 



r+l 



(A + xD) r (B + xC) r 
(x£? + C) r (xA + D) r 

i(r) if(r) 
if(r) i(r) 



= i(l+x) 



+ 2 (1 " X) 



(^-^)r+l = ^ 1+X ) 



Hr) 

-if(r) -i(r) 



+ 2^- X ) 



(A + D) r 
(A + D) r 

if (r) (A-L>) r 

(A — D) 

-if (r) (A - D) r 

(r) ^(r) (A — D) r 

i (r) -if (r) (A + D) r 

if(r) -i(r) (A + D) 



(5.8) 



J(r) 

-iff. 



(5.9) 



The signification of these relations concerning eigenvalues and why we display also (A — D) r+1 
will be explained below. 
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Let us introduce at this point the general possibilities 



i ri+r2 



A r , 



■Br\+r2 A ri 



^4 r2 -|- B ri ® C r2 , 



T21 



D r 

a 



D r , <g> -D 



V1+T2 



= a 



n 



,- 2 "I - Cri 



B, r 



a 



T'Z ■ 



(5.10) 



Since the sequence for odd and even r have some distinct typical features a two-step 
iteration can be of interest for (A ± D). One has (in evident notations) 



2 (A + D) r+2 = (A + £>) 2 ® (A + D) p + (A 
+ (B + C) 2 ®(B + C) r -(B 



D) 2 ®(A-D\ 

C\ ® (5 - C) r 



2 (A - D) 



r+2 



(A + D) 2 ®(A- D) r + (A - D) 2 ®(A + D\ 
(B + C) 2 ®(B- C) r + (B- C) 2 <g) (B + C) r 



(5.11) 
(5.12) 



leading to 

4 (^4 + D) r+2 



(1 + x) 



!(r) 



f (r) 



(r) 



(r) J M 



'(r) 



f(r) 



+ (i-xr 



(r) 



-K 



-l 



(r) 



(r) A (r) 

-K (r) 

-I(r) 
I(r) 



K 



(r) 
f(r) 



(A + D) 





-/, 



(r) 



K, 



(r) 



(r) 



(r) 



(r) 



-if, 



(r) J W 



+ (l + x)(l-x) 



I(r) 





K {r ) 





-K (r) 





-I(r) 








I(r) 





K (r) 










-I{r) 





(A + D\ 



(A + D) r 




(A-D) 









(A + I>) r 






(A + D) 





(A-D) 









(A + D) 




(A + D) r 




(5.13) 







(A + D\ 



(A-D) r 

(A - D\ 



From (5.12) one obtains an analogous result for (A — D) r+2 . 

We now extract one fundamental consequence of the recursion relations (5.8-13). Start- 



1 
1 

(with x-independent matrices X) 



ing with (A + D) 1 = (1 + x) 



and (A - D) 1 = (1 - x) 







these relations imply 



(A + D) r = (1 + x) r X (r , ) + (1 + x) r - 2 (1 - x) z X (r _ 2 , 2) 



+ • • • + (1 + x) z (1 - x) r " 2 X (2 , r _ 2) + (1 - x) r X {0 



for even r and 



(5.14) 

(A + D) r = (1 + x) r X M) + • • • + (1 + x) (1 - x) r - x X (1 , r _ 1) (5.15) 

for odd r, powers of (1 — x) being always even in both cases. Correspondingly (again 
with x-independent 2 r x 2 r matrices Y) 



(A - D) r = (1 + xf- 1 (1 - x) Y (r _ 1;1) + (1 + x)^ (1 - x) a Y (r _ 3 ,3) + • • • + 



r-3 



(i + x) 5 (i- x y 



(S,r-S), 



(5.16) 
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-l) r ). There being here only odd powers of (1 — x). When these series 



where 5 — \ (1 ^ 
are inserted in (5.8) and (5.9) the features above are conserved. 
From the general constraints (see Appendix) 

(A ± D) {A' ± D') = (A' ± D') (A±D), 



(5.17) 



where A = A (x), A' = A (x'), etc. it already follows that the X (resp. Y) must commute. 
Thus with X( r _ 2p ,2p) = X(p), Y( r _ 2 p-i,2 P +i) = Y (p) one must have, for all (p,q) 



[X(p),X ( g)] = [Y(p),Y ( ,)] =0. 
But our recursion relations imply stronger constraints. From (5.8) 



X 



(r+l,2p) - ^ 



l(r) 



K, 



(r) 



K( r ) I(r) 



X (r,2p) + 



f(r) 



(r) 



— K(r) —I{r) 



(r,2p-l) 



Note that 



f(r) 



-AT. 



(r) 



(r) 
f(r) 



= 0. 



(5.18) 



(5.19) 



(5.20) 



Along with recursion relations, systematically exploiting the constraints (A. 2-12) of Ap- 
pendix, one obtains that not only do X( p ), X( g ) commute for each r but 



X( p )X( g ) = X(q)X( p ) = 0, p 7^ g. 



Analogously one can show 



(5.21) 



(5.22) 



These are indeed sufficient and necessary conditions for the eigenvalue spectrum derived 
for r = 1, 2, 3, 4 in sec. 4. 

Let |v p ) denote an eigenstate of X( p ) with the eigenvalue X( p ) |v p ) = v p |v p ) (f p 7^ 0). 
Then for q 7^ p, 

X(,) |v p ) = u^X^X^) |v p ) = 0. 



Hence 

(A + £>) r |v p ) = (1 + x) r - 2p (1 - x) 2p Wp |v p > = 0, p = 0, 2, . . . 
But there is still one more class of constraints. Starting with (see (5.2)) 

Tr(T 1 ) = Tr((A + J D) 1 ) = 2(l+x) 

and noting that (see (5.8)) 

Tr (T r+1 ) = Tr [(A + D) r+1 ) = (1 + x) Tr ((A + D) r ) , 



(5.23) 
(5.24) 

(5.25) 

(5.26) 
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one obtains 

Tr (T r ) = Tr ((A + D) r ) = 2 (1 + x) r . (5.27) 

Hence 

J2v P = 0- (5.28) 

How is this constraint implemented for each r? In our examples (sec. 4, r = 1, 2, 3, 4) we 
saw that 

1. v has a multiplicity 2, saturating (5.27). 

2. f p (p ^ 0) comes with multiplicity, each case providing a subset of zero sum. For 
each p one has, one or more, subsets 

5>W = 0. ( 5 -29) 

i 

3. Here r-th roots of unity play a crucial role. Typically in (5.29) one has 

(1 + x) r ~ 2p (1 - x) 2p e 1 ^* J = 0. (5.30) 

\fc=0 / 

4. For r = 3 one has only cube roots of unity. For r = 4 one has both 2-plets and 
4-plets, (square roots of unity being also fourth roots). 

Let us now consider possible submultiplets from a more general point of view. The 
even and odd subspaces introduced in sec. 4 (even and odd multiplicities of the index 
1 distinguishing them) can be generalized to all r, each one closed under the action of 
(A + D) r and of dimension 2 r ~ 1 . In each one there is exactly 1 state with eigenvalue 
(l + x) r saturating (5.27). Hence one can now consider separately two base spaces of 
dimension (2 r_1 — 1). When r is a prime number (say L) there is a relative simplicity 
concerning the multiplet structure. A Theorem of Fermat (see Ref. 1, Appendix B: 
Encounter with a theorem of Fermat) adapted to our case assures 

2 L - 1 -1 = Z-L, (5.31) 

where / is an integer. Thus for L = 3, 5, 7, 11, etc., / = 1,3,9,93, etc.. Hence an integer 
number of L-plets can span adequately the 2 L ~ 1 — 1 dimensional space with ^I=o e 1_ ^ fc = 0. 
When r is not a prime number each prime factor of r (4 = 2 x 2, 6 = 2 x 3, etc.) can lead 
to submultiplets with zero sum. Finally, if a singlet occurs in {e} the even subspace (i.e. 
apart from (1 + x) r ) it must occur in {o} the odd one with an opposite sign (ex: (1 — x) 4 
in {e} and — (1 — x) 4 in {o} for r = 4). The number of possibilities increase rapidly 
with r. Our study remains incomplete concerning the precise number of multiplets and 
the multiplicity for each higher r. We have however delineated completely, for all r the 
dependence of the eigenvalue spectrum on x (or the spectral parameter 6). 
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Leu us note one point. We changed over from R(0) to R(x) in sec. 3. Thus gives 
conveniently a single parameter x on the anti-diagonal. With the original normalization 



(1 + X) r ~ 2p (1 - X) 2P » e (r^P)m[ + 1 ) e+2pm[-^_ (533) 

Finally, starting with Ti or T 2 (of (4.8)) and implementing recursions one can show 
that the sum of the elements in each row (and each column) of T r is (1 +x) r . The 
sum of basic components of {e} and {o} (|ei), \oi) for r = 3,4 in sec. 4 and their 
direct generalizations) thus each corresponds to the eigenvalue (1 + x) r . The sum of 
these two furnishes the total trace of T r , namely 2 (1 + x) r . All the other states together 
contribute zero trace. Moreover, since for our choice of domains (sec. 3) always < x < 1, 
(1 + x) > 1, (1 — x) < 1 assuming roots of unity in (5.24) for v p . (1 + x) r is the largest 
eigenvalue. This is significant in statistical models. 

In this context one should note that the special status of X^o) in the iterative struc- 
ture. From (5.19) one obtains the matrices structure 



2X( r+ i j0 ) 



X( r ,o) -^(r)X( rj0 ) 
-^<r)X( r ,0) X(r,0) 



(5.33) 



Thus, in each subspace, one can construct iteratively exclusively X( r;0 ) staring from X( 10 ) 
and the corresponding 2 r_1 mutually orthogonal eigenstates: 

1. One with eigenvalue (1 + x) r . 

2. 2 r ~ 1 — 1 with eigenvalue zero. 

The latter provide non-zero eigenvalues for X( p ) with p non zero. One thus obtains the 
complete basis of eigenstates. 

6 Generalizations (n > 2) 

One has for all n 

PPij } = \ {{ji) ® (u) + (f) ® (fj) + e [{ft) ® (fj) + (ji) <g> {ij)}} , 

PP | } = \ ® + ^ + e fa') + ® (u)]} • (6-1) 
These lead (for fundamental blocks or r = 1) to 

T - = (ji) + aj7 ) (ji) , T ?J = (ji) + a£> (ji) , 

T fj = a<7> (ji) + aj } (ji) , T,, = a£> (.//) + a£> (.//) , (6.2) 

where = | ^e mi J 61 ± e m ^ e ^j. For > m^- (resp. mf < ra\ - •*) all Boltzmann 

weights are nonnegative for 6 > (resp. 6 < 0). Note that the fundamental blocks are now 
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(2n x 2n) matrices with only two non-zero elements each. This number does not change 
with n. In a compact notation with indices a e {1, 2, ... , 2n} and a = {2n, 2n — 1, . . . , 1} 
correspondingly, one can write, for r = 1, 



T ab = a { + ] (ba)+a ( J(ba) 



(-) 



with (not only aff = cl^', but) 



,(±) 



a 



(±) 



a 



(±) 



(±) 

6a 



The coproduct rules gives the iterative structure 



^ba 



- ab 



oT (r) 

" ' L cb 



,(+) 



,H f/vi 



)) 



cb ■ 



(6.3) 



(6.4) 



(6.5) 



The fact that, as in (6.3), only diagonal blocks have diagonal elements can be easily shown 
to lead to 



3 r = Tr(T (r) ) =2 [J2 e 



rmL +) 6> 



J=l 



This is the direct multiparametric generalization of the 4x4 case 

3 r = 2e rm ^ e . 



(6.6) 



(6.7) 



It is instructive to study the case n = 2 explicitly. Denoting (as in Ref. 1 with 1 = 2, 
2 = 1) 



°(±) 

h±) 



1(^>.± ,-&>.), C(±) = l(, 



TO 22 )61 -J- g r 



) 



(6.8) 



for r = 1 (with now i — 1,2, i — 4,3 below - a change of notation convenient for displaying 
symmetries) one has 





a + 






















T n = 


















, T22 — 






d+ 





d_ 













a_ 
























a_ 













Tn = 


















j T 2 2 = 








d + 


d_ 






a + 















































Tl 2 = 


c+ 













c_ 


T l2 = 




c+ 






c_ 
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b + 















6_ 



-21 









b_ 
































b + 









(6.9) 



The remaining 8 blocks are given by 

(a+, b+, c+, d+) ^ (a_, 6_, c_, d_) 

One has 



(6.10) 



3f! = Tr (T n + T 22 + T„ + T 2 5) = 2 ((a+ + a_) + (d+ + d_)) = 2 (>" )e + e m - 9 ) . 

(6.11) 

The recursion 



leads to 



% + i = Tr ((a + + a_) (T<? + Tg) + (d + + d_) (t« + 4>)) 



(6.12) 



(6.13) 



a particular case of (6.6). For n = 2 the diagonal blocks, for example, have the iterative 
structures 
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r T (r) 
r T (r) 



a_T^ 

c_T^ 

c_T^ 

a^T { ^ 



•22 



6+TS5 6_T^ 
d + Tg d_T| 



12 



6_T 



12 



(6.14) 



Tjj +1 ^ and T^ +1 ^ are now obtained by setting respectively in T^ +1 ^ and T^ 1 "* a-T~ 
for a + Tjj and so on, systematically in an evident fashion. Their sum gives the transfer 
matrix of order (r + 1) exhibiting the iterative structure. For r = 1 the transfer matrix 
is directly diagonal for all n giving directly the eigenvalues. For n — 2, for example 
(consistently with (6.13)) 



T (i) 



e m<+>« 










































e m<+>* 



(6.15) 



with evident generalization for n > 2. For n = 1 we have systematically explored the 
remarkable structure of the transfer matrix for all r (see (5.14-18)) and consequences for 
eigenstates. A parallel study for n > 1 is beyond the scope of this paper. Our results in 
the section already indicate how the multiparametric aspects start playing on essential 
role. 
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In sec. 5, starting with the 2x2 matrix K (5.3-7) at the level of r = 1 and implement- 
ing tensor products, powerful recursion relations were obtained. We started by relating 
(A, B, C, D) among themselves. For N > 2 one can similarly relate (for a given pair of 
indices the quartet (T^ , T^, T-, T^) given, for r = 1, by (6.2). Evidently one 

can relate through constant matrices only blocks involving the same pair of parameters 

m^j. For this one introduces the matrix 



E ((*) + (*)) 



(6.16) 



generalizing K = (11) + (11) 



1 

1 



Thus for n = 2, one has 



K®K 



1 
10 
10 
10 



The generalization is evident. From (6.8) 



(T n , T„, Tn, T n ) = (T lT , T n , T n , T„) , 
(T n , T n , Tn, T n ) = (T„, T u , T n , T lT ) 



(6.17) 



(6.18) 



with exactly analogous results for other subsets. For each n one has, in evident notations, 
for fixed 



^ ^Tjj, Ty, Tjj, Tjj^) ^Tjj, Tjj, Ty, Ty^) , 
z^W (tv. t^ T - T- l = (T--- T - t- t A 



(6.19) 



We will not attempt to explore in the present paper the applications of such relations 
generalizing our results of sec. 5. 



7 Spin chain Hamiltonians 

Our construction of odd dimensional Hamiltonians (sec. 4, Ref. 1) can be adapted to the 
present even dimensional cases as follows. One has, taking derivatives and setting 9 = 0, 
for n = 1, 











X 







X 








X 


x+ 





X 








x+ 
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where x-t 



setting (starting from (6.8)) 



.(+) 



(-) 



and R (0) is obtained by setting 9 = in -4i? (9). For n = 2 



6± 



1 

2 
1 

2 



l 12 



l i2' ) y • 



from (2.7) and 



R(0) 



where 

Dn = Dii 
An = in 



a_, 6_, 6_, a_ 



d ± 
c± ■ 



±m { 22 ] 



hi 



hi 











An 





D 22 


A22 








A-22 


D-22 





An 











diag. 
anti-diag. 



D22 = D22 = [c + , d + , d + , c + 
A22 = A22 = (c_,d_,d_ 



(7.2) 



(7.3) 



diag. 

c 



• ,.(7-4) 
anti-diag. 



The extension of our formalism for n > 2 is straightforward. 

For r sites the standard result for the Hamiltonian is (see sources cited in Ref. 1^ 



R(0) 



k,k+l 



(7.5) 



k=l 



where for circular boundary conditions fc + l = r + l«l. We intend to present a more 
complete study of our spin chain elsewhere. But here one already sees how the two aspects, 

multistate (higher spins at each site) and multiparameter (jn\^^j get directly associated 

for our hierarchy. In our previous papers [HIE] and here again (see (3.4), (3.5)) we showed 
how the passage to imaginary parameters can lead to unitary R(9). The corresponding 

R (0) has only an overall factor i which can be extracted from the sum (7.5). 



8 Remarks 

I. Status of eigenstates for n = 1: For the simplest 4x4 braid matrix in our 
hierarchy the construction of the eigenvalues and eigenfuntions of transfer matrices of 
successive orders (r = 1, 2, 3, 4, etc.) has attained the following stage: 

(1) The transfer matrix at the level r has been expressed in the form of a 2 r x 2 r matrix 

T,= £ (l + xr 2p (l-xf p X (p) , (8.1) 

p=0,l,2,...,p m 



17 



where x = tanh| (^n^ ~ m ii^ 9 an d 2p m = r (resp. r — 1) for r even (resp. 

odd) (2p m = r — (1 — ( — l) r ) /2). being the spectral parameter, m 1 ^ two free 
parameters (see (3.3)). The matrices X( p ) are constant ones (x-independent) and 
satisfy 

X(p)X ( ,) = X ((?) X (p) = 0, p ^ g. (8.2) 
They can be computed systematically via the recursion relations (for r — >■ r + 1). 



(2) For 



X(p) |p w > = u^i) |p (i) > , (u^i) ^ 0) , 

X (<?) \P(i)) = °> P ^ ? ( 8 - 3 ) 



and f denote phase factors which come in multiplets of zero sum formed by roots 
of unity corresponding to r and its prime factors (see examples in sec. 4). The index 
"i" denotes such possible multiplicity of each p. The exception of the zero sum rule 
corresponds to p — 0. One obtains for each r, twice i>( ) = 1 giving 

Tr (T r ) = 2(1+ x) r (8.4) 

a general constraint obtained via recursions. This multiplicity 2 corresponds to two 
2 r_1 dimensional subspaces (see "even", "odd" subspaces defined in sec. 4) each 
providing just one eigenstate of X( ) with non-zero eigenvalue (1 + x) r . 

(3) Thus the problem has been reduced to construction of eigenstates of each X( p ) sep- 

arately, reducing the dimension by considering each (even, odd) subspace by turn. 
This involves solving sets of linear constraints with only positive and negative inte- 
gers as coefficients. One finally keeps only the non-zero eigenvalues for each p , they 
being associated with zero eigenvalues for X( 9 ), q ^ p. In fact as noted below (5.33), 
it suffices to construct the full sect of mutually orthogonal eigenstates of X( rj0 ), all 
but one in each subspace having eigenvalue zero. 

(4) As already stated (sec. 5) our results remain incomplete concerning the pattern of 

possible multiplets and submultiplets corresponding to roots of unity provided by r 
and its prime factors - and multiplicities of such multiplets. A canonical enumeration 
when r has a very large number of prime factors seems to be an unlikely possibility. 
Nor have we established rigorously that in (8.3) are always ±1 or higher roots of 
unity phases factors. This what happens in examples (r < 4) of sec. 4 and directly 
leads to the following obligatory constraint (8.4). 

(4) We have completely, and for all r, extracted the ^-dependence of eigenvalues in (8.1). 
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II. Comparisons with standard six vertex and eight vertex models: From (2.5) 
and (3.1) our Yang-Baxter matrix (for n = 1) with corresponding normalizations, is 



R(0) = PR (9) 



a + 








a_ 





a_ 


a + 








a + 


a_ 





a_ 








a + 



■5) 



where a± 



g'"-n c 



(-)* 



±e m ii 



and, equivalently, 



1 x 
0x10 
1x0 
x 1 



(8.6) 



where x 



tanh | 



711 



(+) 



11 



— m 



(-) 



n y 6*. Let us now compare this to the very well known 4x4 

six vertex and eight vertex models - concerning which it is sufficient to a cite a standard 
text book [7] and review articles [HI [8] which cite basic sources. All such cases are of the 
form 

a d 
b c 
c 6 
d a 



R{9) 



7) 



For (8.6) a = c = 1, 6 = d = x and for (8.5) a = c = a+, b = d = a_. For six 
vertex, crucially, d = and (a, 6, c) is being given, according to the regime, by circular 
or hyperbolic functions. For eight vertex, famously, elliptic functions appear and d ^ 0. 
For our case d ^ but arguably, one has maximal simplicity and symmetry compatible 
with non-trivial solution for a 4 x 4 Yang-Baxter (or braid) matrix. We have all eight 
vertices but (with, say (8.5)) ( see fig. 1) corresponding to a = c = a + and (see fig. 2) 




fig. 1 
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corresponding to b = d = a_. In six vertex the last two vertices are excluded (d = 0). 
We will not study, in this paper, the implications of the results below (8.7) concerning 
various properties of our model (compare the relevant detailed study of eight vertex in 
Ref. 7). But we would like to contrast our approach to the construction of eigenstates and 
extraction of eigenvalues of T r with that via Bethe ansatz in standard six vertex models 
[5]. The systematic study of RTT constraints are particularly relevant (see Appendix A). 



In six vertex the Bethe ansatz construction involves pushing (A [9] 
the product B (9i) B (62) . . . B (9 r ) acting on one single state 



<8> 



D (9)) through 
1 N 



<8> 



1 11 . . . l) r and eliminating unwanted terms to obtain a complete set of eigenstates corre- 
sponding to the sets of resulting constraints. This involves solving nonlinear equations. 
Our resorts to programs and numerical studies for higher r's. 

In our case recursion relations for (A ± D) r (for r — > r + 1) are sufficient to attain 
the stage systematically presented in part (I) of this section. One solves , at each stage, 
linear equations with integer coefficients. In fact since each r (see (5.5-7)) (B ± C) T = 

1 



K {r) (A ± D] 



, where Km in the 2 r 



x 2 r matrix with 2 r 1 -times K = ^ ^ on diagonal 

the actions of (B ± C) r follows immediately from those of (A ± D) r . They need hardly 
be studied separately. However that Km and hence (B ± C) r connect the even and odd 
subspaces. On the other hand (A + D) cannot be pushed through a product of B (#)'s 
displayed above. The nearest approaches are typically our (A. 10) and (A. 12) (changes 
of signs in (A ± D) and (B ± C) are to be noted as (A + D) is pushed through). Subtle 
analytic properties, unlike six vertex and particularly eight vertex models, play no role in 
our case. One has only to look at the x (or 9) dependence in (8.1). 

In a simple situation the contrast between our model and the standard eight vertex 
one shows up very clearly: For the asymptotic case 9 tending to infinity one sets x = 1 
in (8.6) to obtain 

10 1 
110 
110 
10 1 



R = R 



Whereas, choosing adequately the normalizing factor (see sec. 7 of [9] and in particular 
eq. (7.13)) for the standard case the corresponding limiting form is 



R= (l,g- 



1) 



diag. 



±R. 



•9) 



When diagonalized (or block diagonalized) our classes of matrices, in general, lose braid 
(or Yang-Baxter) property unless the diagonalizer has a corrected tensor structure (sec. 
4 and Addendum of ref. 3). 



III. passage to higher dimensions n > 2: As already emphasized before, a major 
interest of our simple model for n = 1 (the 4x4 braid matrix) is that it is the first one in 
a hierarchy of (2n) 2 x (2n) 2 braid matrices with 2n 2 free parameters at each level. Some 
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of the simplicity of the n — 1 case is inevitably lost as n increases. But we have pointed 
out in sections 6 and 7 how certain basic features vary in a simple, canonical fashion. 
Thus for example, 

(1) 

Tr (TW) = 2 e rm " +>f? j , for all n. (8.10) 

(2) Each blocks of the T matrix has just two non-zero elements (see (6.2) and (6.9)) 

out of (2n x 2n), namely (9) for all n. 

(3) For a fixed pair of indices the blocks (TV,-, Tjj, T^) can be quite simply 

related among themselves (see (6.18)) via direct generalizations of the matrix K = 

^ q for the 4x4 case. For n — 1 such relations led to recursion relations yielding 
(8.1-3). 

(4) Spin chain Hamiltonians present (see eqs. (7.1-5)) a simple canonical sequence as n 

increases. 

We hope to study the higher dimensional cases more fully elsewhere. 

Acknowledgments: One of us (BA) wants to thank Pierre Collet and Paul Sorba for 
precious help. 
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Appendix A 

RTT constraints (n = 1) 



From (3.10) one obtains with K from (3.9) 



A B 
C D 



= x 



A' B' 
C D' 

B' A' 

d' a 

B A 
D C 



A' B' 
C D' 



B A 

D C 

B' A' 

jy a 



A B 
C D 

C D 
A B 

a d' 

A' B' 



a d' 

A' B' 

C D 
A B 



(A.l) 



The last step follows from x ^ x' on both sides of the first two expressions since under 
this interchange (see (3.7)) x" — > — x". The consistency of the last two implies that each 
element of the total matrix on the right (apart from the factor x") must be symmetric 
in (x, x'). This can indeed be verified starting from r = 1,2,3,... using the standard 
construction in sec. 3. 

From the last two steps (the factor x" cancelling) one obtains, with only upper or 
lower signs, 

(A ± D) (A' ± D') = (A' ± D') (A±D), 
{B ± C) (B' ± C) = (B' ± C) (B±C), 
(A ± D) (B' ± C) = (A' ± D') (B±C), 
(B ± C) {A' ± D') = (B' ± C) (A ± D) . (A.2) 

One gets 8 relations of the type 

(.1 / ; .\/', - M[M 2 ) = x" (M^M 4 - M 5 Mq) = x" (M 3 M' 4 - M^Mq) , (A.3) 

with 

(Mi, M 2 , M 3 , M 4 , M 5 , M 6 ) 

= {(A, B, B, A, C, D) , (A, C, B, D, C, A) , (B, A, A, B, D, C) (B, D, A, C, D, B) , 
(C, A, D, B, A, C) , (C, D, D, C, A, B) , (D, B, C, A, B, D) , (D, C, C, D, B, A)} (A.4) 



In sec. 5 (see (5.4-7)) we obtained, for any x and all r, 

(B, C) = K {r) (A, D) , (A, D) = K {r) (5, C) , 



(A.5) 



where (x-independent ) matrix is given. Hence multiplying in (A.3) by on the 
left one gets another set with (Mi, M 3 , M 5 ) replaced by (Mi, M 3 , M 5 ), where thus 



(A, B, B, A, C, D) — ► (B, B, A, A, D, D) 



(A.6) 
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and so on. For another class of relations we introduce 

f<^I(x» ± l) (A.7) 

From (3.7) 

f « = coth n (9 - 9') , f (_) = -cosech/i (9 - 9') (A.8) 
One can show that (with upper or lower signs) 

(A ± D) (A 1 t D') = f W (B ± C) (B 1 t C) + f (_) (B' ± C) (B T C) , (A.9) 
(A ± D) (B' t C) = f (+) (B ± C) {A' t D') + fH (B' ± C) (A^D). (A.10) 

Again, as for (A. 5), (A. 6), multiplying the above from the left by one gets another 
set of relations such that 

{B ± C) (B' ± C) = f (+) {A ± D) {A' t D') + fH (A' ± D') (At D) . (A.ll) 

and so on. For "two-steep" relations redefine x" with indices Xj ? - = ^ Xj and f/. + ^ = 

coth/x (6*j — f^ - " 1 =-cosech/i (9; L — 9j) correspondingly. One obtains, for example, for 
arguments (xi,x 2 ,x 3 ) 

(A + D\ (B - C) 2 (B + C% = f x { 2 +) 4 +) (5 + C), (B - C) 2 (A + D) 3 + 

fi?& (S + C) 1 (5-C) 3 (A + J D) 2 + 
f^fW (5 + C) 2 (5 - C), {A + L>) 3 + 
ff-^ (B + C) 2 (S - C% (A + £>)! . (A.12) 

See the relevant remarks in sec. 8. 
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